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Abstract 

Within an effective field theory framework, we obtain an expression for the next-to- 
leading term in the 1/m expansion of the singlet QQ QCD potential in terms of Wilson 
loops, which holds beyond perturbation theory. The ambiguities in the definition of 
the QCD potential beyond leading order in 1/m are discussed and a specific expression 
for the 1/m potential is given. We explicitly evaluate this expression at one loop and 
compare the outcome with the existing perturbative results. On general grounds we 
show that for quenched QED and fully Abelian-like models this expression exactly 
vanishes. 
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1 Introduction 



After the discovery of the first heavy-quark bound states, the ip and the T systems, it was 
soon realized that a non-relativistic picture seemed to hold for them. This is characterized by, 
at least, three scales: hard (the mass m, of the heavy quarks), soft (the relative momentum 
of the heavy-quark-antiquark |p| ~ mv, t) < 1), and ultrasoft (US, the typical binding 
energy E ~ mv 2 of the bound state system). It was also seen that, if one wanted to describe 
the whole spectrum of the if) and the T systems, a perturbative evaluation of the potential 
was not sufficient. This triggered the investigation of these systems by all sorts of potential 
models (see [IJ for some reviews), which are, in general, quite successful phenomenologically. 
Since then, a lot of effort has been devoted to obtaining the relevant potentials to be used 
in the Schrodinger equation of such models from QCD [Q, |3|, [|, || ^| [7| |8|, ^| by relating these 
potentials with some Wilson loops that could eventually be computed on the lattice or by 
using some vacuum models. The expression for the leading spin-independent potential, of 
O(l/m ), has been known since long and corresponds to the static Wilson loop @, [|. The 
expressions for the leading spin- dependent potentials in the 1/m expansion, of 0(l/m 2 ), 
have been calculated in Refs. @, ||. The 1/m corrections to these potentials have proved to 
be very difficult to obtain. To our knowledge the spin-independent case has been addressed 
only in Refs. || (|. The result of Ref. |J does not reproduce the one-loop perturbative QCD 
potential (see discussion at the end of section 5.1) and, therefore, appears to be incomplete. 
In Ref. the author does not succeed in obtaining suitable finite expressions. We conclude, 
therefore, that the question of the 1/m corrections to the QCD potential has not been 
settled yet and hence deserves further studies. In this work we will present an ab initio and 
systematic calculation of the QCD potential up to 0{l/m). We will get a new expression of 
the 1/m potential that is finite, consistent with one-loop perturbative QCD and suitable to 
be evaluated by lattice simulations. 

We will perform the calculation by integrating out in two steps the hard and the soft scales 
characterizing the heavy-quark-antiquark system. This is implemented by introducing suit- 
able effective field theories. This approach allows us to express the heavy-quark-antiquark 
dynamics in terms of systematic and controlled expansions. It has proved to be a powerful 
computational tool in several different situations. For instance, the hard log corrections 
(~ lnm) to the Eichten-Feinberg-Gromes potentials were computed in this way in Ref. 
|T0| (see also [§]). Moreover, we believe that the effective field theories provide a suitable 
framework where eventually some long-standing conceptual questions will be clarified. In 
particular, the extent of validity of the naive potential picture for the heavy quarkonium dy- 
namics, assumed in potential models, could be affected by the consideration of extra degrees 
of freedom such as hybrids and pions. 

The two QCD effective field theories that arise from integrating out the scales m and 
mv are called NRQCD and pNRQCD respectively. Non-relativistic QCD (NRQCD) was 
first introduced in Ref. [IT]. It has an ultraviolet cut-off much smaller than the mass m and 
much larger than any other scale (in particular much larger than Aqcd, which means that the 
matching from QCD to NRQCD can always be done perturbatively [|l|, |13|). NRQCD has 
proved to be extremely successful in studying QQ systems near threshold. The Lagrangian of 
NRQCD is organized in powers of 1/m, making in this way explicit the non-relativistic nature 
of the described systems. The maximum size of each term can be estimated by assigning 
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the soft scale to any dimensionful object. In order to connect NRQCD with a potential 
picture the degrees of freedom with energies much larger than mv 2 have to be integrated 
out. Once this is done, one is left with a new QCD effective field theory called potential 
NRQCD (pNRQCD) 0, [TJ[]. Strictly speaking, pNRQCD has two ultraviolet cut-offs, Ai 
and A2. The former fulfils the relation mv 2 <C Ai <C mv and is the cut-off of the energy of 
the quarks, and of the energy and the momentum of the gluons, whereas the latter fulfils 
mv <A 2 Cm and is the cut-off of the relative momentum of the quark-ant iquark system, 
p. This theory has been thoroughly studied, and its matching with NRQCD performed, in 
the situation mv 3> Aqcd- in this case the matching can be performed perturbatively. In this 
paper we will allow Aqcd to be as large as mv and, therefore, we cannot rely on perturbation 
theory. Nevertheless, we will assume that the matching between NRQCD and pNRQCD can 
be performed order by order in the 1/m expansion. We will present, for the general situation 
Aqcd <j m v, the matching of NRQCD to pNRQCD at the next-to-leading order in the 1/m 
expansion in the singlet sector (to be defined later). This will prove to be equivalent to 
computing the 1/m potential. Formulas that are similar to those in the classical papers of 
Refs. [||, H |6| will be worked out. No further degrees of freedom with US energy besides the 
singlet will be considered. This means that non-potential effects will be neglected in this 
paper (in the perturbative situation this would be equivalent to working at zero order in the 
multipole expansion). A detailed study of the matching between NRQCD and pNRQCD in 
the situation Aqcd < ttiv, including US effects, will be worked out elsewhere. 

The paper is organized in the following way. In section ^| we introduce NRQCD up to 
order 1/m and discuss its static limit. Moreover, we define what will be pNRQCD in the 
present context. In sections ^ and [| we derive the 1/m corrections to the potential, by 
matching NRQCD to pNRQCD. The Green functions are worked out in the Wilson loop 
language. The 1/m potential can be written in a simple way as insertions of chromoelectric 
fields on a static Wilson loop. In section || we compute the potential perturbatively up to 
one loop. Quenched QED and Gaussian models of the QCD long-range dynamics are also 
discussed. Finally, section 6 is devoted to the conclusions and the Appendix to show how 
unitary transformations affect the form of the potential. 



2 NRQCD and pNRQCD 

The Lagrangian of NRQCD up to order 1/m reads 

-\g%g^\ (1) 

where ip is the Pauli spinor field that annihilates the fermion, \ is the Pauli spinor field that 
creates the antifermion, iD° = id® — gA°, and zD = zV + gA. The matching coefficients 
Cp ~ 1 + 0(a s ) are not going to be relevant here. For simplicity, light fermions are not 
explicitly displayed. Their inclusion does not change the results of this paper, although it 
changes the expression of some intermediate formulas. 

1 For related work on these issues within the perturbative regime of NRQCD we refer to [[i~2| |l6], |l7| . 



2 



The Hamiltonian associated to the Lagrangian ([!]) is 



mi wi2 

H (0) = J rf 3 x I (gaga + gaga) ^ 

if™ = -i / (d 2 + £ C «<r ■ B) ^ = \ J (d 2 + gcfcr ■ B) X , 

and the physical states are constrained to satisfy the Gauss law. We are interested in the 
one-quark-one-antiquark sector of the Fock space. 

It will prove convenient to study the static limit. The one-quark-one-antiquark sector of 
the Fock space can be spanned by 

|n;x 1; x 2 )( ) := ^ t (x 1 )x(x 2 )|n;x 1 ,x 2 )(°), V Xl ,x 2 (2) 

where |n;x 1 ,x 2 )(°) is a gauge-invariant eigenstate (up to a phase), as a consequence of the 
Gauss law, of with energy £^°)(xi,x 2 ); |n;xi,x 2 )w encodes the gluonic content of the 
state, namely it is annihilated by x^( x ) an d VK X ) (^ x )- ft transforms as a 3 Xl ® 3* 2 under 
colour SU(3). The normalizations are taken as follows 

(0) (m;xi,x 2 |n;xi,x 2 ) (0) = 5 nm (3) 

(°)(m; x 1; x 2 |n; y 1; y 2 >(°) = <W (3) (*i - yi)5 (3) (x 2 - y 8 ) ■ (4) 

Notice that since does not contain any fermion field, |n; x 1; x 2 )(°) itself is also an eigen- 
state of if (0) with energy E^°)(xi,x 2 ). We have made it explicit that the positions Xi and 
x 2 of the quark and antiquark respectively are good quantum numbers for the static solu- 
tion |n; x 1; x 2 )(°\ whereas n generically denotes the remaining quantum numbers, which are 
classified by the irreducible representations of the symmetry group (substituting the 
parity generator by CP). The ground-state energy £q (x 1; x 2 ) is usually associated to the 
static potential (see [JTSJ for important qualifications on this association) and the remaining 
energies £^°)(xi, x 2 ), n / 0, are usually called gluonic excitations between static quarks. 



They can be computed on the lattice (see for instance |L8|). Translational invariance implies 
that E^'(xi, x 2 ) = E^\r), where r = Xi — x 2 . 

The gap between different states at fixed r will depend on the dimensionless parameter 
Aq CD r. In a general situation, there will be a set of states {n us } such that E^' (r) ~ mv 2 
for the typical r of the actual physical system. We denote these states as ultrasoft. The 
aim of pNRQCD is to describe the behaviour of the ultrasoft states (pNRQCD has been 
introduced in |TJ| and discussed in different kinematic situations in [fOJ). Therefore, all 
the physical degrees of freedom with energies larger than mv 2 will be integrated out from 
NRQCD in order to obtain pNRQCD. In the perturbative situation Aqcd 7 " 1, which 
has been studied in detail in fll5fl , {n us } corresponds to a heavy-quark-antiquark state, in 
either a singlet or an octet configuration, plus gluons and light fermions, all of them with 
energies of 0(mv 2 ). In a non-perturbative situation, which we will generically denote by 
Aqcd?" ~ 1, it is not so clear what {n ns } is. One can think of different possibilities. Each 
of them will give, in principle, different predictions and, therefore, it should be possible to 
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experimentally discriminate between them. In particular, one could consider the situation 
where, because of a mass gap in QCD, the energy splitting between the ground state and 
the first gluonic excitation is larger than mv 2 , and, because of chiral symmetry breaking of 
QCD, Goldstone bosons (pions/kaons) appear. Hence, in this situation, {n us } would be the 
ultrasoft excitations about the static ground state (i.e. the solutions of the corresponding 
Schrodinger equation), which will be named the singlet, plus the Goldstone bosons. If one 
switches off the light fermions (pure gluodynamics), only the singlet survives and pNRQCD 
becomes totally equivalent to a quantum-mechanical Hamiltonian, thus providing us with 
a qualitative explanation of how potential models emerge from QCD. In addition, we shall 
see below how quantitative formulae can be provided in order to calculate the potentials in 
QCD. 

In this paper, we will only study the pure singlet sector with no reference to further 
ultrasoft degrees of freedom. In this situation, pNRQCD only describes the ultrasoft excita- 
tions about the static ground state of NRQCD. In terms of static NRQCD eigenstates, this 
means that only |0; Xj, x 2 )^ is kept as an explicit degree of freedom whereas |n; x l5 x 2 )(°) 
with n ^ are integrated out0. This provides the only dynamical degree of freedom of the 
theory. It is described by means of a bilinear colour singlet field, S , (x 1 ,x 2 ,i), which has 
the same quantum numbers and transformation properties under symmetries as the static 
ground state of NRQCD in the one-quark-one-antiquark sector. In the above situation, the 
Lagrangian of pNRQCD reads 

sUid - /i s (x 1 ,x 2 ,pi,p 2 ) )S, (5) 



^pNRQCD 

where h s is the Hamiltonian of the singlet (in fact, h s is only a function of r, p 1; p 2 , which is 
analytic in the two last operators but typically contains non-analyticities in r), pi = — iV Xl 
and p 2 = — ?V X2 . It has the following expansion up to order 1/m 

h s (x u x 2 , Pl , pa) = J£- + p- + V (r) + (— + —) V x {r). (6) 
Zrrii 2m2 \mi m 2 / 

In this work we will present the matching between NRQCD and pNRQCD within an 
expansion in 1/m using the static limit solution as a starting point. Whereas this can be 
justified within a perturbative framework, in the non-perturbative case the validity of the 
1/m expansion cannot be generally guaranteed^. We note that this assumption is implicit 
in all the attempts at deriving the non-perturbative potentials from QCD we are aware 
of. Furthermore, we would like to emphasize the following important point. The matching 
calculation can be done independently of what the specific counting in pNRQCD is. As 



argued in [13|, [L7|, [15[ , when doing a matching calculation we are integrating out high energy 
degrees of freedom. Hence, the Wilson coefficients (potentials) that this integration produces 
are independent of what the low energy dynamics is. The 1/m expansion just provides a 
convenient way to organize the matching calculation. One should not conclude that the 
relative size of the different potentials in the pNRQCD Lagrangian is trivially dictated by 
the 1/m expansion. 



2 In fact, we are only integrating out states with energies larger than mv 2 and all the states with n^O 
will be understood in this way throughout the paper. 

3 See p9| for an example where certain degrees of freedom cannot be integrated out in the 1/m expansion. 
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3 The Wilson loop matching 

In this section we will work out the matching up to order 1/m, in a language that is close 
to the traditional approach of Eichten and Feinberg ||, [|. The matching between NRQCD 
and pNRQCD is done by enforcing suitable Green functions to be equal in the two theories 
at the desired order in 1/m. In this section we shall consider space-time Green functions, 
which are more conventional in non-perturbative studies. 

Let us consider an interpolating state with a non- vanishing overlap with the ground state, 

^ t (xi)0(xi,x 2 )x(x 2 )|vac), (7) 

where <fi may in principle be everything that makes the above state overlap with the ground 
state \0)^ . We will use here the following popular choice^: 

0(y,x;t) =Pexp ji0jTtfo(y-x) ■ A(x - s(x - y), t)| , (8) 

where P is a path-ordering operator. We also define 0(y,x;t = 0) = 0(y,x). It is going to 
be useful to introduce a n (xi,x 2 ), defined by 

^ t (x 1 )0(x 1 , x 2 )x(x 2 ) |vac) = fl n(xi, x 2 ) |n; x 1; x 2 ) (0) . 

n 

The states |n)( ) being defined up to a phase, it is convenient to fix them in such a way that all 
the coefficients a n are real. The identification of the singlet from the state (0) depends on the 
interpolating field used in NRQCD. This dependence is reflected in different normalization 
factors Z. The matching condition reads 

X t (x 2 ,t)0(x 2 ,xi;t)^(x 1 ,t) = Z 1/2 (xi,x 2 ,pi,p 2 )5(x 1 ,x 2 ,t). (9) 

As in h s , the normalization factor Z only depends on r, pi and p 2 , and is given in the form 
of an expansion in 1/m as 

Z(xi, x 2 , Pl , pa) = Z (r) +(— + —) Z x {r) + iZ ltP {r)r . ( J* _ -El) + . . . , (10) 

\mi m 2 / Vmi m 2 / 

where we have made use of the fact that the NRQCD Lagrangian (|IJ) as well as the pNRQCD 
Lagrangian @ are invariant under a CP plus m x «-> m 2 transformation. 
We will match the Green function Gnrqcd defined by 

Gnrqcd = (vac| X t (x 2 ,T/2)0(x 2 ,x 1 ;T/2)^(x 1 ,T/2) 

x^(yi, -T/2)0(y 1; y 2 ; -T/2) X (y 2 , -T/2)|vac), (11) 

with the corresponding Green function in pNRQCD (here and in the rest of the paper, if not 
explicitly stated, dependence on x ii2 , p ii2 is understood) 

Gpnrqcd = Z^ e~ iTh ° ZtV^(3)( Xl _ yi)5 (3)( X2 _ y2 ). (12) 

4 With this choice we assume that the ground state has the £+ quantum numbers (see |p|). This is so in 
perturbative QCD as well as non-perturbative QCD according to the available lattice simulations. 
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3.1 Calculation in NRQCD 

We expand Gnrqcd order by order in 1/m 

r - m 1 r^ (1 ' 0) a. 1 _i_ 

<^NRQCD — ^ NRQCD ~>~ _ ^ NRQCD ~t~ ^NRQCD ~~r • • • • 

Tfl\ iit'2 

Integrating out the fermion fields one gets 

GfiacD = (^□)5 (3) (xi-yi)5 (3) (x 2 -y 2 ), (13) 
G{!r°qcd = y^ 2 / f(^t)} n 6^\^ 1 -y 1 )5^^ 2 -y 2 ). (14) 

Analogous formulas hold here and in the following for G^rqcd- For simplicity we will not 
display them. The angular bracket (...) stands for the average value over the Yang-Mills 
action, Wa is the rectangular static Wilson loop 



W a = Pexp jf tfeM^z) j 



the angular bracket (...)□ = (. . . Wq) stands for the average over the gauge fields in the 
presence of the static Wilson loop ((l)n = (Wa))- For further convenience we also define 
((...))□ = (. . . W n ) I (Wa). We use the convention that, if not specified, fields act on the first 
quark line, e.g. D(t) = D(xi,£), E(t) = E(xi,t), and so on. We have used time reversal: 
(B(i))a = — (B(— t))a to eliminate the spin-dependent term in Eq. fllip . Therefore, we can 
already state that no spin- dependent potential appears at 0(l/m) M. 

It is useful to introduce at this point two identities involving covariant derivatives and 
Schwinger lines: 

i) D(x, tU(t, x, t', x) = (pit, x, t', x)D(x, t') + ig tdt" <j>(t, x, t", x)E(x, t"U(t", x, t', x), 

Jt' 

ii) D(xi,t)0(xi,x 2 ;t) = igr x / ds s </>(xi, x'(s); t)B(x'(s), £)0(x'(s), x 2 ; t) 

Jo 

+0(xi,x 2 ;t)V Xl , x'(s) = sxj. + (1 - s)x 2 , 

where 0(t,x, t',x) = Pexp |— ig J^dt" Aq(x, t")j- In the first Ref. ||], both identities i) 

and ii) were derived. Identity ii) corrects their Eq. (4.7c). As a by-product of the above 
identities, we get 

rT/2 



V xl (W D > = i0 / ^(E(t)) D + (O i (-T/2)) D -(0 / (T/2)) n , (15) 

J-T/2 



where 



Oi(t) = igr x / ds s 0(x 1; x'(s); f)B(x'(s), i) </>(x'(s), x 2 ; £), x'(s) = x 2 + sr, 
Jo 

0/(t) = igrx f ds (1 - s) 0(x 2 , x"(s); t)B(x"(s), t) 0(x"(s), x x ; t), x"(s) = xi - sr. 

w 



,T/2 




Li? 1 









Let us note that time-reversal symmetry gives {ig~E(t)) n = (igHj(—t))n and (Oj(— T/2)) n = 
-<0,(T/2)) D . 

Using the above relations, Eq. ([TJ]) can be worked out giving 

GfeSU = ^^ 1 (^) + ^(%)v^+r(o / (T/2).o j (-r/2)) D 

/(T/2) ■ E(f)) D - foj^jfc (| + t) (E(t) • 0,(-T/2)) D 
+!L dt ^|t-t'|(E(t)-E(t'))nU (3) (xi-yi)5 (3) (x2-y 2 ). (16) 

2 J-T/2 J 

3.2 Calculation in pNRQCD 

Let us consider the Green function defined in Eq. flT2|) . Expanding it up to order 1/m, we 
obtain 

rr _ ^(0) , 1 ^(1,0) 1 ^(0,1) 

'j'pNRQCD — ^pNRQCD "+* ~~ ^pNRQCD "+* ~~ ^pNRQCD- 

nh\ r/t-2 
Inserting Eqs. @ and (0) into Eq. (0), we obtain 

cSrqcd = Z e-^ 3 )(x 1 -y 1 )5 (3) (x 2 -y 2 ) (17) 
^pNrqcd - z oe s — ~ ~ I v r 



\ Z 2 V ' Z 

l T (VV )-r^ + zT^- i TV 1 
Z Zjq z 



iT / (VZ ) _ , o (V 2 Z ) (VZ > 



8 1 Zn Zq -^0 



+^ f 2(W ) • v Xl + • (w ) + (Wo) 



4 l-v z 







-?(Wo)4^xi-yi)* (3) ( 



6 v u/ J 



x 2 - y 2/ 



where V = V r . As in the NRQCD case, we do not display the analogous formulas for 
GipNRQCD- In order to keep Eqs. (PT7|) and ( |TSj ) simpler, we have already used the fact that 

Z and Zx, p can be chosen as real. This follows from the matching to G^rqcd dH) once any 
constant phase in Z is conventionally set to zero. 

3.3 Matching 

At O(l/m ) we match Eq. (P) with Eq. (0). We obtain 

Vq = lim ^\n(W n ). (19) 
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At 0(l/m) we match Eq. flTjjp with Eq. (|TB|). We obtain 



{ i 1 ~ |) « /( T / 2 ) • E (*)»° " (* + |) « E W ■ 4 (-T/2)» D } 



4 J-T/2 

-\((O f {T 12)0^/2)%-^- f T/2 dt f T/2 dt'\t-t'\((E(t).E(t'))}X (20) 

2 4i J-T/2 J-T/2 / 



A similar expression with neither end-point string contributions nor normalization factors 



has been obtained in |J. The right-hand side of Eqs. ( |T9| ) and ( |20|) can be shown to be 
finite in the T — > oo limit. This is not obvious from the point of view of a pure Wilson-loop 
calculation in the spirit of Ref. j§], as is apparent from the difficulties met by the author of 
H (the normalization factors are crucial in order to get a finite expression). For this reason, 
and because such kind of arguments may become relevant to future analyses of Wilson-loop 
correlators, we mention here the relevant steps of the proof. 

a) Inserting the identity operator J2 \n)^ ^(n\ into the Wilson- loop average, the latter may 
be written as 



<wb> = xy 



iE { ° ] T „2 



a 



n ' 



Doing the same for the average (E(t)) n one obtains 

(E(t)) D = E e "^° >T °n {0) (n\E\n)^ + £ e-^+^^-^t^ ®(n\E\m)®. 

b) End-point strings containing the operators Of and Oj select intermediate states with 
quantum numbers different from the singlet. This can be checked directly on the definitions 
of these operators. As an immediate consequence, correlators containing the operators Of 
and Oi in Eq. (p0|) vanish in the T — > oo limit and do not contribute to the potential. 

c) From Eq. ( |ToD and time-inversion invariance of the chromoelectric field, it follows that 

-(Wo) = (0) ((%E|0>(°\ 

(°)((%E|n)W 



[VZ ) = 2J2 °0«n 



E n — Eg 



d) Inserting the identity operator into the correlator (E(t) ■ E(i')) n , it may be written as 
(E(t) ■ E(0)d = E Mm (0) (n|E| S )(°) (o) (s |E|m)( )e-^ 0) +^ 0) >^ 



n,m,s 



for £ > With the above points a)-d) it is easy to show that the right-hand side of Eqs. 
(|T9|) and ( |2"0"D is finite in the large-T limit. Their explicit expressions read 

V = (21) 
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9_ 
AT 




+ i 



T {VZ 
4 Z n 



rT/2 


r T/2 n 


/ dt 


/ dt'\ 


l-T/2 J 


-T/2 



dt'|t-t'|((E(t)-E(t')»n 

2 

(°)(n|^E|0)(°) " 



E, 



(o) 



p(o) 

-C/n 



(0) 



n| 5 E|0) 



(o) 



i 0) ) 2 ' 



(22) 



Analogously, we can obtain matching expressions for the normalization factors. At 0(l/m°), 
the result is Z Q = \a \ 2 . At 0(l/m) a more complicated equality is obtained, which involves 
a combination of Z\ and The fact that we have two equations and three independent 
functions: Z 1; Z lp and V\ at 0(l/m) makes the determination of V± ambiguous. This 
ambiguity is intrinsic in the sense that any value of Z\ iP and V± that fulfils Eq. ( p2[) will lead 
to the same physics (as far as one consistently works at higher orders in 1/m), and has to 
do with the fact that a quantum-mechanical Hamiltonian is defined up to time-independent 
unitary transformations (see Appendix). As a consequence, the ambiguity in the definition 
of V\ also gives rise to ambiguities in the definition of the potentials of order 0(l/m 2 ) and 
higher. In the next section we will fix Z\ iP (and then V\) by imposing an extra matching 
condition, which will prove to be particularly convenient. 



4 The quantum-mechanical matching 

In the previous section we have done the matching comparing Green functions in NRQCD and 
in pNRQCD. In this section the comparison is between states and matrix elements in NRQCD 
and in pNRQCD. The calculation, in terms of states, will be closer in philosophy to the usual 
quantum mechanics calculations in perturbation theory |20| (see also [21]). Moreover, the 
whole procedure will share some similarities to the adiabatic approximation and the Born- 



Oppenheimer method as used in atomic physics calculations |22| . The underlying assumption 
is that the difference of energies among states labelled with different n is much larger than 
the difference of energies among states labelled with the same n. In our case, since we only 
aim at correctly reproducing the ground state spectrum, we only need that the splitting 
between the ground state and the first gluonic excitation be larger than the typical splitting 
of the states of the ground state (taken of 0{mv 2 ) by definition). This is nothing but the 
condition we have assumed throughout the present paper. 

H is not diagonal in the basis of (\n)^) with respect to the n labelling. We consider 
instead a basis of states, labelled as 

|n;xi,x 2 ), (23) 
such that the Hamiltonian H is diagonal with respect to these states 

(m;xi,x 2 |-£/ln;yi,y 2 ) = 8 nm E n {yL X , x 2 , pi, p 2 )5 (3) (xi - yi)<5 (3) (x 2 - y 2 ) , (24) 

where E n (x.i, x 2 , pi, p 2 ) is an analytic function in pi, p 2 , and such that they are normalized 

as 

(m;xi,x 2 |n;yi,y 2 ) = 6 nm 5 {3 \x 1 - yi)5 (3) (x 2 - y 2 ) . (25) 
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Conditions fl24|) and 



give 



ff|n;yi,y 2 ) = J rf 3 Xid 3 x 2 |n;x 1 ,x 2 ) J E n (xi,x 2 ,pi,p 2 )5 (3) (x 1 - yi)5 (3) (x 2 - y 2 ) . (26) 

A set of states |n) and an operator E n that satisfy Eqs. fl25|) and (p6|) can be obtained from 
the static solutions In)*- -* and E^ to any desired order of accuracy by working out formulas 
analogous to the ones used in standard quantum mechanics perturbation theory []2D| . If we 
write |n) and E n as an expansion in 1/m, 



|n) = |n)(°) + — In)^ ) + — (n)^ 1 ) + . . 
mi m 2 

E n = E^ + J_£(i.o) + J_ s (o,i) + ... 
mi m 2 



(27) 



(28) 



we obtain at 0(l/m) for n = (when not specified, states and energies are calculated in Xi 
and x 2 ): 



10} 



(1,0) 



E, 



(0) 



(°)(n|£/E|0)( ) > 



1 



p(o) 



E 



2 V' X1 E^-E^ 



1 W(n[gg|j)W W(j|^E|0)W . (0) 

2 ^ fp(o) p(»)\/p(») P (o)x ' 



(29) 



^^(X! - y a )* W (*a - y 2 ) = (°)(0|^|0> 



(0) 



V X1 ' 



E 

n^0 



(0) 



<n|(/E|0) 



(0) 



P (0) p(0) 



5( 3 )(xi- yi )5( 3 )(x 2 -y 2 ). (30) 



These equations may be derived from the identities 

a) (0) <n|D x »(°) = V xl , 
which follows from symmetry considerations and 

b) w= w yy VB5tJ , 

which follows from explicit calculation. Analogous formulas hold for the antiparticle contri- 
bution. 
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4.1 Matching 

The aim of pNRQCD is to describe the behaviour of |0). The integration of high excitations 
is trivial using the basis ( p3|) since, in this case, they are decoupled from |0). Therefore, the 
matching of NRQCD to pNRQCD is basically to rename things in a way such that pNRQCD 
reproduces the matrix elements of NRQCD for the ground state. The matching conditions 
in this formalism read as follows: 

|0) = S^vac) and S (xi, x 2 , pi, p 2 ) = h s (x.i, x 2 , pi, p 2 ) . (31) 

Using Eq. (|3T| ) and Eq. (|9|) we obtain the normalization factor, given the interpolating field, 

^(xx, x 2 , Pl , p 2 )5 (3) (x! - yi)5 (3) (x 2 - y 2 ) = (vac| X t (x 2 )0(x 2 , x 1 )^(x 1 )|0; y 1; y 2 ). (32) 

It is explicit now that h s and, hence, the potential, fixed by Eq. (|3l|) , depend neither on 
the normalization factors nor on the specific shape of the end-point strings. Moreover, Eqs. 



(|3T|) and (|32|) provide matching conditions at any finite order in 1/m. 

From Eqs. fl31~D and (^) one could draw the conclusion that the ambiguity in the com- 
putation of Vi, discussed in the previous section, has disappeared. This is not the case. 
Actually, Eqs. Q2"7| ) and ( PBQ with Eqs. ( ]2T)D and (^) only give one of the possible solutions 
of Eqs. (^) and (|25|). Indeed, these equations do not completely fix the state |n). In stan- 
dard quantum mechanics the state is fixed up to an arbitrary constant phase. In our case, 
since we only diagonalize in the n space, this phase becomes a unitary operator^ given a 



state |n) that fulfils Eqs. ( p26|) and fl25|), this means that the states 

J d 3 x.[ d 3 x' 2 |n, x£, ^ 2 )e l °^^'^ ) 5^\^ 1 - Xl )5 (3) (x 2 - x 2 ), 
with = O n , still fulfil Eq. ( ^5l) and the Hamiltonian is still diagonal in n with 



E n ^e iu -E n e 



iO n 771 — iO„ 



This freedom reflects the fact that h s is defined by the matching up to an arbitrary unitary 
transformation or, which is the same, up to an arbitrary unitary field redefinition: 

and so does the Z. In order to make calculations easier, we have taken advantage of this 
freedom by fixing the relative phase between |0) and lO)^ -* following the standard choice of 
quantum mechanics perturbation theory. With this choice, we obtain Eqs. (^) and (|30|); as 
we will see below, this will allow us to obtain a compact expression for V\ in terms of Wilson 
loops. Furthermore, the procedure can be easily generalized at any finite order in 1/m. 

5 Note that the phase is not completely arbitrary if we demand the state |n) to coincide with the un- 
perturbed state |n)(°) in the limit 1/m — > 0. This constrains O n (x,p) to smoothly go to zero in the limit 
1/m -> 0. 
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From Eqs. fl29|) and (^) we can perform the matching at 0(l/m). We obtain 



7 i/2 

7 l/2 



do, 



E 



(33) 



^0 -C'O 



(0) 



p(o) 

-Cm 



v xi - 



(0) 



(n|(/E|0) 



(or 



(0) 



(n| 5 E|0) 



(0)' 



4 0) 



(4 0) - S 



(Oh 2 



:v xiJ e 



(0)> 



(°)(n|^E|j}( ) ( )(j|^E|0)(°) 



,(o) Wz? (o) 



v 1 / 2 



E / ,-,(0) , .IU! , .|UJ , ,11 

C°)(n|^E|0}(°) 



(oh 



2E a 

n^0 



(0) 



From Eqs. and ( |3"1~D we get, up to 0(l/m) 

Vn(r) - ^ (0) < 



Vl(r) 



E 

n^0 



(0) 



H<?E|0> 



(o) 



(0) 



E, 



(o) 



(34) 
(35) 

(36) 
(37) 



At first sight, this result seems of limited practical utility, since it depends on the exact 
and complete solution of the bound state at order l/m°. Fortunately this is not the case. In 
fact, Eq. ( pOD suggests a very simple form for the 1/m potential, 



Vi 



T^oo 



( 4T„ 


,772 


fT/2 


/ dt 






'-T/2 J 


-T/2 



(E(f) • E(t')»o - «*)»□ • «E(f)»c 



(38) 



is 



Using the same technique as outlined in aj-<i) of Sec. 3.3, it can be proved that Eq. (|38| 
finite and equal to Eq. fl3~7[) . 

Equation (|38|) is the main result of the present work. It gives the leading 1/m correc- 
tion to the heavy-quark potential expressed in terms of chromoelectric field insertions in a 
static Wilson loop, thus avoiding the explicit computation of the normalization factors. The 
potential appears, therefore, in the same form as the potential calculated in Refs. 0, |^, || 
and is suitable for lattice evaluations similar to those performed in ||. Moreover, the above 
expression expanded in a s gives the full perturbative series of the 1/m potential in the regime 
where this expansion makes sense (i.e. mv ^> Aqcd)- I n the next section we will calculate 
from it the leading non- vanishing perturbative contribution to the 1/m potential, which will 
coincide with the results of Refs. [^|, ^4| . For the calculation of some higher order pertur- 
bative corrections to the V\ potential we refer to [p5| . Finally, it is interesting to note that 
Eq. (^) also holds if the correlators are evaluated on a Wilson loop with infinitely large 
time extension (i.e. if we make the limit T — > 00 in the Wilson loops for t and t' finite, and 
then evaluate the integrals). This observation greatly simplifies the perturbative calculation 
and may also do so for non-perturbative ones. 



5 Applications 

In this section we will consider Eq. (|38D in perturbative QCD and in quenched QED. In 
the first case we will evaluate the contribution to the potential at the leading non-vanishing 
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order in a s . In the second we will show that the potential V\ vanishes exactly. This may 
be relevant to several Gaussian models used in the phenomenology of the non-perturbative 
dynamics of the strong interaction. We will shortly comment on this. 



a) 




b) 




Figure 1: Dashed lines indicate Coulombic exchanges, curly lines transverse gluons. Graph 
a) is the disconnected graph cancelling the T 2 (VVo) 2 /12 term in Eq. Graph b) (and 

the symmetric graph which is understood) gives the only non-vanishing contribution to the 
1/m potential at order a 2 in the Coulomb gauge. In other gauges additional diagrams may 
be present. 



5.1 Perturbative QCD 

We perform the calculation in the Coulomb gauge. Since Eq. fl3"8p is a gauge-independent 
quantity, the result will hold in any gauge. At order a s , the right-hand side of Eq. (|38|) only 
gives a self-energy type of contribution, since the chromoelectric fields are inserted on the 
same quark line. The first non-vanishing contribution to the potential (i.e. depending on r) 
appears at order a 2 . Three types of diagrams arise: 

i) Disconnected graphs. An example is shown by graph a) of Fig. [I]. These graphs cancel in 
the difference of the right-hand side of Eq. (^) . 

ii) Connected graphs with end-point strings contributions. Non-disconnected graphs involving 
gluons attached to the end-point strings vanish in Eq. (|38[) . As noticed in the previous section 
the T — > oo limit can also be performed in the correlators before doing the time integrals in 
Eq. fl3~8|). This turns out to be quite useful here, making the irrelevance of these graphs to 
the potential manifest. 

Hi) Connected graphs with no end-point strings contributions. The only non- vanishing graph 
contributing to Eq. fl38D is graph b) of Fig. [l]. The analogous graph involving the triple- 
gluon vertex, but with two transverse gluons attached to the chromoelectric fields on the 
quark line, can be shown to vanish in the limit of Eq. (|38|) by explicit calculation. 
Therefore, considering all relevant contributions (which reduce to graph b) of Fig. [I]) we get 
from Eq. ( p8|) at order a 2 

Vr" = ~C F C A ^. (39) 
Equation (|39|) coincides with the result of Refs. J23 , |24j (which may also be obtained using 



the rules of and with the non-Abelian part of the perturbative calculation of the 1/m 
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potential done in [25]. In Ref. [25 also an Abelian contribution to the 1/m potential is 
found. We have stressed throughout the paper that there is no unique way to define the 
1/m potential and, therefore, different matching procedures are, in general, expected to give 



different expressions for it. One may suspect that this is due to the fact that in Refs. p3[ ^4 



the Coulomb gauge is used whereas in [£5J the Feynman gauge is used. However, this has 
more to do with the freedom we have in quantum mechanics to change the form of the 
Hamiltonian by a unitary transformation without changing the physics, rather than with 
the gauge fixing dependence itself. What typically happens in perturbative calculations, 
which match non gauge invariant Green functions, is that depending on the gauge one uses, 
one gets a different form of the potential. The different forms are equivalent and can be 
obtained from each other by unitary transformations. More precisely, from the discussion in 
the Appendix it follows that V^ pe /m, at O(a^), can be rewritten in terms of an 0(a s /m 2 ) 
potential However, whereas the Abelian piece found in [26 can be obtained from the 



general expressions for the 0(l/m 2 ) potential ||, this is not so for the non- Abelian piece in 
Eq. (|39|), since there is not non- Abelian contribution at tree level. Therefore, we conclude 
that Eq. (|38|) represents a genuine new potential not considered in the past in the non- 
perturbative parametrizations of the QCD potential in terms of Wilson loops. Finally, let us 
mention that, to our knowledge, Vi eTt has never been computed using Wilson loops before. 

5.2 Quenched QED 

Let us consider quenched QED. In this situation the Wilson-loop average is exactly known 
(it reduces to Gaussian integrals, see for instance P7[| ) and can be written as 



f Q 2 f i f i fT /2 r T/2 

{W n ) = exp { -2- / ds ds' dt dt'nrjfijit-t', 

{ 2 Jo JO J -T/2 J -T/2 



[s - s')r) , (40) 



where fij(t,r) = fji(t, r) = (Ei(t, r)Ej(0, 0)}. As a direct consequence (for a derivation, see 
for instance ||28|| ) we obtain 

«E(t) ■ E(f )»□ - «E(t)» D • «E(f)»D = f*(t - f , 0). (41) 

Therefore, since the term fu(t — t',0) is a self-energy-type contribution, the potential con- 
tribution of V\ (as defined in Eq. ( ]38|) ) vanishes exactly in quenched QED. The same result 
is obtained in Ref. RTTf . 



Several QCD vacuum models |], 0, [29| seem to approximate the Wilson-loop long-range 
non-perturbative dynamics with expressions analogous to Eq. (flOl). The above considerations 



made for quenched QED may be relevant to them. As a matter of fact, 1/m corrections do not 
seem to show up there. Therefore, it is tempting to generalize the exact result of quenched 
QED to Gaussian models of the QCD vacuum. Nevertheless some words of caution are 
needed. A model consists in an approximation of QCD, which is supposed to coincide with a 
relevant limit on the QCD dynamics. This limit is unknown by definition. Therefore, how to 
implement it is, in the practice of several models, not well established. While it is clear that 

6 For the general, unequal mass, case the Abelian term is of the type l/(m 1 + m 2 ) and the field redefinition 
discussed in the Appendix would transform it into a l/(m\m-2) term. 
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starting from a Gaussian expression for the Wilson loop a relation like Eq. ( f4T|) is going to 
hold, it is not guaranteed that 1/m potentials will not show up if the Gaussian approximation 
of the Wilson loop is implemented at some intermediate step of the calculation]]. 



6 Conclusions 

We have obtained an exact and non-perturbative expression for V\, the QQ potential at order 
1/m. This expression is shown in Eq. fl38l). The perturbative calculation of it up to order 
a 2 agrees with [^, |24| and is consistent with the one-loop contribution of [2?J. However, 
the existence of a non-perturbative contribution at order 1/m was, to our knowledge, not 
considered before in the literature, be it in phenomenological applications |IJ or in attempts 
at obtaining the non-perturbative potentials from QCD |4], |5], [J. We note that via a unitary 
transformation the 1/m terms can, in certain circumstances, be reshuffled in 1/m 2 (and 
higher) terms in the potential (see Appendix). When these circumstances apply (i.e. when 
V\ <C m 2 v 2 ), our results do not have immediate consequences on phenomenological models 
where the full set of 1/m 2 potentials are considered and input as ansatze functions. However, 
when they do not (i.e. when V\ ~ m 2 v 2 ), our results imply that a 1/m potential should be 
included in those models. In either case, they are extremely important for the attempts at 
obtaining the potentials from lattice QCD, since, upon doing the above-mentioned reshuf- 
fling, the 1/m potential given by Eq. ( |38"D generates 1/m 2 terms different from the ones 
calculated so far ||. 

The most promising way to have an estimate of the non-perturbative behaviour of the 
correlators appearing in Eq. fl38|) is by a lattice simulation. In practice this can already be 
done from the lattice data of Ref. ||, since the correlators we get are of the same type. Such 
lattice data could be of interest for at least two reasons. If the potential found on the lattice 
happens to be small in the long range, it would support the Abelian dominance picture. In 
the short range it should show the interplay of the perturbative (known) region with the 
non-perturbative one. 

When there are no additional US degrees of freedom, the adopted non-perturbative pro- 
cedure can be easily generalized to the evaluation of higher-order terms in the inverse mass 
expansion JHJ) . The inclusion of US degrees of freedom in the non-perturbative regime will 
be discussed elsewhere. 
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7 The reader may, for instance, compare the formulas of Refs. || with the predictions of some Gaussian 
models discussed in Ref. (291. 
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Appendix A 



In this appendix we show that there exists a unitary transformation that reshuffles momen- 
tum-independent 1/m terms in 1/m 2 momentum-dependent and momentum-independent 
terms. Let us consider the Hamiltonian 



H = f + V,,(r) + m. (42) 
2m m 



The unitary transformation 



U = exp( {W(r),p} , 



m 



transforms H — > H' = W H U . More explicitly, under the condition {W, p}Cm (which is 
necessary in order to maintain the standard form of the leading terms in the Hamiltonian, 
i.e. a kinetic term plus a velocity independent potential) H' reads 

h' = ^- + v + — + -w-(w ) + -^w-(Wi) 

2m m m m z 

+^'(W(vy„)) - JLtf, (v«r)}} + o . (43) 



w = -iy r Wn 



2 vvy n 



By choosing 



all the 1/m potentials in the Hamiltonian ( pE3"D disappear. The price to pay is the appearance 
of new terms at order 1/m 2 (and higher). Of course, the leading size of these new terms 
is the same as the original V\jm. In particular, since p, V, ~ mv and Vq ~ mv 2 , we 
get W ~ Vi/(m 2 v 3 ) and {p l , {p 7 , (V i H / - 7 )}}/m 2 ~ V\/m. The size of the remaining 1/m 2 
induced terms is V 2 /(m 3 v 2 ) and hence, it depends on the size of V\/m, which is a priori 
unknown. On general grounds the maximum size of V\ is given by the largest available scale, 
namely mv, and hence at most V\ ~ m 2 v 2 . Reasoning in the same way, the 1/m 2 potentials 



(calculated via the quantum-mechanical matching in [ 30 1 ) are not bigger than mv . However, 
from the condition {W, p} m, it follows that the reshuffling of the 1/m potential to 1/m 2 
potentials may be done in the way above only if V\ m 2 v 2 . As a consequence, all the terms 
of 0(V 2 /(m 3 v 2 )) have a size much smaller than mt> 2 . More specifically, if, for instance, 
V\ ~ m 2 v 3 , then the terms of 0(V 2 /(m 3 v 2 )) are of order mv 4 and hence suppressed by a 
factor v with respect to the original 1/m potential (as well as with respect to the 1/m 2 
potentials obtained from the quantum-mechanical matching) . We note that in perturbation 
theory (v ~ a B ) V\ ~ m 2 v 4 and, therefore the terms of 0(V 2 /(m 3 v 2 )) are suppressed by a 
factor v 2 . 
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